LOOP-ERASED WALKS AND TOTAL POSITIVITY 



SERGEY FOMIN 



Abstract. We consider matrices whose elements enumerate weights of walks 
in planar directed weighted graphs (not necessarily acyclic). These matrices 
are totally nonnegative; more precisely, all their minors are formal power series 
in edge weights with nonnegative coefficients. A combinatorial explanation 
of this phenomenon involves loop-erased walks. Applications include total 
positivity of hitting matrices of Brownian motion in planar domains. 



1. Introduction 

Then [Pooh] dropped two in at once, and leant over the bridge to see which of 
them would come out first; and one of them did; but as they were both the same 
size, he didn't know if it was the one which he wanted to win, or the other one. 

A. A. Milne Chapter VI] 

Consider a stationary Markov process whose state space is a connected planar 
domain f2, or a discrete subset thereof. In the continuous version, assume that 
all trajectories of this process are continuous; in the discrete version, the possible 
transitions from one state to another are described by a planar directed graph. 

Let B be an absorbing subset of the boundary of 17 (i.e., once the process reaches 
a state in B, it stays there), and suppose that the points ai, a2 ^ B and 61, 62 G i? 
are such that any possible trajectory that goes from oi to 62 must intersect any 
trajectory going from 02 to bi . See Figure |^. 
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Figure 1. Some possible locations of ai, 02, 61, 62 



Let us now consider two independent realizations of our Markov process, starting 
at points ai and 02, respectively. Suppose that we know that the other endpoints 
of these trajectories are bi and 62, although we do not know which of the two ended 
up where; neither do we know how much time it took for each one to reach its 
destination. We then ask the usual Bayesian question: which of the two possibilities 
(i.e., ai 61, a2 ^ 62 or ai ^2, 02 ^ 61) is more likely, and by how much? 
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In particular, is it necessarily true that the first match-up is more likely than the 
second one? 

This paper answers the last question — and various related ones — in the affir- 
mative. To make it more precise, let Xij denote the hitting probability (or the 
corresponding density) that the trajectory originating at will hit the target B 
at location bj . Then our claim amounts to the inequality xiia;22 > 2:122^21 • In 

other words, the determinant of the submatrix ( ) of the hitting matrix 

\X21 X22 ) 

X of our Markov process is nonnegative. In fact, the determinant of any square 
submatrix of X can be shown to be nonnegative, so the hitting matrix is totally non- 
negative. To explain this phenomenon, we provide a combinatorial interpretation 
of such determinants in the spirit of the celebrated Karlin-McGregor nonintersect- 
ing path approach; note that this approach cannot be applied directly to the case 
under consideration, since the trajectories are allowed to self-intersect. The crucial 
ingredient of our combinatorial construction is Lawler's concept of loop-erased walk. 

To give a flavor of the main result of this paper, here is its simplest version, 
stated loosely: we prove that the 2x2 determinant xiia;22 — 2:12X21 is equal to 
the probability (or density) that the trajectories that started out at ai and 02 will 
end up at 61 and 62, respectively, and furthermore the first trajectory will never 
intersect the loop-erased part of the second one. 

We give conditions under which the minors of hitting matrices are positive (as 
opposed to nonnegative), so the matrices themselves are totally positive. This 
latter property immediately implies a number of others, for example, simplicity 
and positivity of the spectrum, and the variation-diminishing property. 

The paper is organized as follows. Sections ||-|| are devoted to preliminaries of 
various kinds. Section |^ introduces walk matrices and hitting matrices of directed 
networks. Section previews classical results by Karlin-McGregor and Lindstrom on 
total positivity of walk matrices of acyclic directed networks or associated Markov 
chains. Section ^ gives an account of some of the results obtained in |ll| on 
resistor networks and their Dirichlet-to-Neumann maps; although these results are 
not used in the rest of the paper, they provided the primary motivation for this 
investigation (cf. Acknowledgments below). Section || introduces loop-erased walks. 

The main results are presented in Sections We give combinatorial formulas 
for the minors of walk matrices and hitting matrices of directed networks (see 



Theorems 6T, OA, 7T, and TJl) and Markov processes (see Theorems 1_A and 7.5), 
and apply them to various examples of discrete and continuous Markov processes, 
such as one-dimensional Bernoulli process and two-dimensional Brownian motion. 
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2. Directed networks and associated matrices 

Most of the material in this section is basic graph theory (see, e.g., Sec. 4.7]), 
although some of the terminology is not standard. 

A directed network T = (V, E, w) is a directed graph (loops and multiple edges 
are allowed) with vertex set V and edge set E, together with a family of formal 
variables {w{e)}eeE, the weights of the edges. 

The notation a b will mean "edge e goes from vertex a to vertex b." More 
generally, 

/r) -1 ^ ei 62 es e„, , 

(Z.L ) a = flo — > fli — > a2 — > • • • — > ttm ~ o 

will denote that the edges ei, . . . , form a walk of length m from a to 6 that 



successively passes through oi, 02, . . . . We will sometimes shorten (2.1) to a — > 6, 
where tt denotes the walk. The weight of tt is by definition given by w{t:) = Y[ w{ei). 
The degenerate walk a — > a of length has weight 1. 

No finiteness conditions are imposed on F; we only ask that the number of walks 
of any fixed length between any two vertices a, 6 G ^ is at most countable, so that 
the following formal power series is well defined: 

(2.2) T¥(a,6) =^^«;(ei)---u;(e™), 

m 

where the second summation goes over all walks of length m from a to 6, as given 



in (2.1). In other words, W{a,b) is simply the generating function for the weights 
of all walks from a to b. The matrix W — {W{a, b))a^bev will be called the walk 
matrix of the network T. (The term "Green function" would perhaps fit better, 
but it is already overused.) 

One easily sees that W = (I — Q)^^, where 

Q ^ {Q{a,b))a,bev , Q{a,b) ^^w{e) , 

a — >b 

is the "weighted adjacency matrix" of the network, and / is the identity map/matrix. 
The identity W = QW + I shows that for any b eV, the function Wb : a W{a, b) 
is "Q-harmonic" on F \ {6}, i.e., satisfies Wb{a) ~ Q(a, c)Wb{c), for a^b. 

Example 2.1. See Figure 0. 
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Figure 2. Directed network and its walk matrix 



In most applications, the weight variables w{e) are specialized to numerical (usu- 
ally nonnegative real) values. One should then be careful while introducing the walk 
matrix W ^ since the power series ( |2.2| ) may easily diverge. The following example 
is typical in this regard. 
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Example 2.2. Consider a Markov chain with a countable set of states V. When- 
ever the transition probability p{a, b) from state a to state b is nonzero, let a 6 
be an edge of E, and let its weight be w(e) — t ■ p{a, b), where t is a formal vari- 
able. Then W = Em=o*'"-F"" = (I - tP)~\ where P = (p(a, 6))a,bev is the 
transition matrix of the Markov chain. Note that we cannot simply put t = 1 and 
w(e) = p{a, b) for a b, since this may lead to divergence of the power series 
for W{a,b). If, however, this power series does converge, then 

expected number of times the process 
visits b provided it starts at a 



(2.3) W{a,b) 



Example 2.3 (Bernoulli random walk). Let V = Z, with edges connecting n and 
n + 1 (in both directions) for all n E Z. Let the weights be given by 



(2.4) w{e) 



\p if n — > n -|- 1; 
[q if n + 1 n, 

where p + q = 1 and p > q > 0. It is classically known (see, e.g., Section 1, 
(11)]) that 

'(p-g)-i iffc>0; 

k 



(2.5) W{n,n + k) 



(p-g)-i(^) iffc<0. 



Hitting matrices. In what follows, we will also need a variation of the notion of a 
walk matrix, which arises in situations where T comes equipped with a distinguished 
subset of vertices dT C V, called boundary of T. Accordingly, intT = V — dT is the 
interior of F. 

For a G V and b £ dT, let X{a, b) be the weight generating function for the 
walks from a to 6 of nonzero length whose internal vertices all lie in the interior 
of F. Formally, 



m>0 

where the second sum is over all walks ( |2.1| ) such that ^ 9F for < i < m. 
The hitting matrix X is then defined hy X ^ {X{a,b))aev,bedT ■ Also, the term 
"hitting matrix" will sometimes be used for submatrices of X, with the notation 
Xa,b = {X{a, b))aeA,b^B tor Ac V, Be dV. 

Example 2.4. See Figure 




A = {01,02} 
5 = {61, 62} 



Xa,b — 



9495 91939497 



1 



1 - 919293 V 91^29596 919697 / 



5F = {6i,62} 



Figure 3. Hitting matrix of a directed network 
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Example 2.5. Let us continue with Example p.2| . Unlike in the case of walk 
matrices, it is now permissible to set t — 1. Then X(a,b) is nothing but the hitting 
probability; more precisely, it is the probability that the Markov process with initial 
state a visits the boundary after leaving a, and furthermore the first boundary state 
it hits is b. 

The matrix XgT,dr can be expressed in terms of the weighted adjacency matrix Q 
using the notion of Schur complement. Recall that for a matrix M of block structure 

M = ( ) , the Schur complement M/F is defined by M/F = C - DF'^E. 



^E F ^ 

The following statement is a straightforward corollary of the definitions. 
Proposition 2.6. We have 

(2.6) / - Xar^av = (/ - Q)I{I - Q^nt^ant^) , 

where Qintr,intr = iQia,b))a,beintr ■ (By a common abuse of notation, in formula 
(p^), letter / denotes identity matrices of three different sizes.) 

3. Acyclic graphs and total positivity 

This section reviews some classical results concerning directed networks whose 
underlying graph is acyclic. Although these results will not be used in subsequent 
proofs, they will serve as inspiration in the study of the general (non-acyclic) case. 

For an acyclic network F, let the boundary dV be the set of sinks in F: 

dr^{aeV:^a^b}. 

Under these assumptions, the hitting matrix A is a submatrix of the walk matrix 
of F, and the entries A(a, b) of A enumerate the paths from each vertex a to each 
sink b with respect to their weight. 

For an arbitrary pair of totally ordered subsets A = {ai,...,ak} C V and 
B — {bi, . . . C dT of equal cardinality fc, consider the corresponding k x k 
minor of the hitting matrix: 



(3.1) det(A^,B) = det A(a„ 6,),=i,...,fc 

V 3 = 1. -.k 

The following fundamental observation goes back to Karlin and McGregor ||l^ (in 
the case of Markov chains) and Lindstrom ||l^ . 

Theorem 3.1. In an acyclic network, the minors of the hitting matrix are given by 

(3.2) det(AA,B) = ^ sgn((T) ^ w{Tri) ■ ■ ■ w{TTk) , 

'^^^1' non- intersecting path families 

oil * ''ct(I)' ■ ■ ■ ; <^k > 

where the first sum is over all permutations a in the symmetric group Sk (inter- 
preted as bisections A —> B), and the second sum runs over all families of vertex- 
disjoint paths which connect the vertices in A to the sinks in B , assigned to them 
by a. (Recall that the weight w{'k) of a path is by definition equal to the product of 
the weights of its edges.) 



Theorem B.l, which is proved by a path- switching argument, has many applica- 
tions in the theory of stochastic processes, enumerative combinatorics, and beyond 
(see, e.g., ||l| and references therein). Most of these applications are based on the 
following corollary, treating the case of planar acyclic networks. 
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Corollary 3.2. Let be a planar simply- eonnected domain whose boundary dfl 
is a simple closed Jordan curve. Let T be an acyclic network, embedded into fl in 
such a way that the edges of T do not intersect. Suppose that A = {ai, . . . ,afc} 
and B — {61, . . . , b^} C dV lie on the boundary dQ, and furthermore these vertices 
appear in the order ai, . . . , a^, 6^, . . . , 61 if traced counter-clockwise. Then 

(3.3) det{XA,B) = 

non-intersecting patli families 
ai — Oi, ... , afc — > Ofe 

In particular, if the weights of all edges are nonnegative, then dci{X a,b) > 0, and 
furthermore Aet{X a' ,b') > for any subsets A' C A and B' C B. 

Thus any such matrix Xa,b is totally nonnegative (cf., e.g., j^, 10 ). 

4. Resistor networks and Dirichlet-to-Neumann maps (after [TT]) 

This section reviews some well known and some fairly recent results from dis- 
crete potential theory — or, equivalently, the theory of resistor networks, — following 
D. V. Ingerman's thesis ijlll and the paper ||] by E. Curtis, D. V. Ingerman, and 
J. Morrow. 

A resistor network is essentially an undirected graph (V,-E) (i.e., a graph in 

which each edge a 6 is paired with an edge b ^ a) together with a conductivity 
function 7 satisfying 7(e) = 7(e')- The conductivities 7(e) can be viewed as cither 
positive reals, or as independent variables taking positive real values. 
All our networks will be presumed connected. 

For a "potential" function m : ^ R, the corresponding "current" function 
/„ (more precisely, the function that gives currents out of each vertex, for a given 
collection of potentials) is given by 

/u(a) = 5]7(e)(w(a)-w(6)). 

a — 'b 

The Kirchhoff matrix K = {K{a, b))a,bev of a resistor network represents the linear 
map u ^ lu from potentials on V to the corresponding currents. Thus 



K{a,b) 



a — *h 



^ 7(e) if a = 6. 

e , 

— >c,c^a 

Suppose that the set of vertices V is partitioned into two disjoint subsets dT and 
intV , as above. The potential functions u : — > M satisfying Kirchhoff 's Law 

I^{a) =YK{a,b)u{b) = , for a e intT , 

b 

are called ^-harmonic. In other words, the value of u at an interior node a should 
equal the weighted average of its values at the neighbors of a. 

The values of a 7-harmonic function u at the interior nodes of T are uniquely 
determined by the values of u at the boundary nodes. This allows us to define 
the Dirichlet-to-Neumann map A of the resistor network as the map that sends a 
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function / : dT ^ K to the current out of the boundary nodes of the unique 7- 
harmonic continuation of /. The Dirichlet-to-Neumann map is represented by the 
response matrix A = (A(a, 6))a,be9r of the network. A straightforward calculation 
yields the following formula. 

Proposition 4.1. Theorem 3.2] The response matrix A is the Schur comple- 
ment in the Kirchhoff matrix K : 

A = K/ KintT,intT , 

where Ki„tr,intr is the submatrix {K{a, b))aMeintr ■ 

For further discussion of Dirichlet-to-Neumann maps, see jH] and references 
therein. 



The comparison of Propositions 4.1 and |2.6| shows that the relationship between 
the Kirchhoff matrix and the response matrix is very similar to the relationship be- 
tween the weighted adjacency matrix and the hitting matrix XgY.dr- This analogy 
is not accidental: we will see very soon that the response matrix of a resistor net- 
work is closely related to the hitting matrix of the Markov chain associated to the 
network in the standard way |^. Recall that this Markov chain has the vertices 
in V as its states, and the numbers 

p(a,6)- 



as transition probabilities. The following statement is a straightforward conse- 
quence of the definitions. 

Proposition 4.2. (6.6)] The hitting matrix Xgr^ar of the Markov chain asso- 
ciated to a resistor network is related to network 's response matrix A by 

(4.1) Xar,ar = 1- Kq^A, 

where Kq — {KQ{a,b))a.bedr is the diagonal part of a principal submatrix of the 
Kirchhoff matrix: Kq^o, b) — K{a, b)5ab- 



Formula (4.1) shows that the hitting matrix of a resistor network is, up to renor- 
malization, its response matrix, that is, the matrix of the Dirichlet-to-Neumann 
map. 

Next, let us look at the minors of the response (or hitting) matrix of a resis- 
tor network. Since the underlying directed graph of the network is almost never 

acyclic (indeed, any two adjacent vertices a and b give rise to a cycle a b ^ a), 
the nonintersecting-path formulas from Section ^ do not apply. Instead, we will 
refer to a remarkable determinantal formula discovered by D. Ingerman (see [|[ 
Lemma 4.1]), reproduced below without proof. 

Let A = {oi, . . . , flfe} and B — {bi, . . . ,bk} he two disjoint ordered subsets of dT 
of the same cardinality k. Let 

Aa,s = (A(ai,6j))i=i,...,fc 
j=i,...,k 

denote the corresponding submatrix of the response matrix A. 
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Theorem 4.3. (D. Ingerman) The minor det{AA,B) of A is given by 
(4.2) 

where, for a permutation a € Sk , we denote 
(4.3) 

CAA,B) = L={n,,...,TTk) : P^ths ai ^ . . . , ^ 6.(,) are 1 

[ vertex-disjoint and stay in the interior of F J 

and K-jj denotes the submatrix of K whose rows and columns are labelled by the 
interior vertices that do not lie on any of the paths in n. 

For a connected resistor network with positive conductivities, the submatrix 
KintT,intT IS always invertible (cf. Proposition 4.1), so the denominator in (4.2) 
does not vanish. 

It follows from Proposition |4.2| that the corresponding minors of the hitting 
matrix of the network are given by 

det(X^,i3) = „ -det(A^,B) 

(4.4) 

E w{'Ki) ■ ■ ■ w{TTk) dei{K^) 
sgnlfij 2^ Aei{K,ntT^ntT)\\r,aAK{a,a)' 

under the assumptions of Theorem 4.3. 

Just as in Section || — and for exactly the same reasons, — formulas (4^) and (4.4) 
simplify considerably in the case of a planar resistor network, yielding the following 
results. 

Corollary 4.4. Assume that a connected resistor network with positive conductiv- 
ities is embedded into a simply- connected planar domain fl in such a way that its 
edges do not intersect and the boundary vertices in dV are located on the topological 
boundary of Q. Let A = {ai, . . . , Ofe} and B — {bi, . . . ,bk} be two disjoint subsets 
of boundary vertices of the same cardinality k; we presume that the vertices in A 
are ordered clockwise, while those in B are ordered counter-clockwise. Then 

(4.5) det(A^,B) = (-1)''^ ^(^i)---^(^O- , w!?^^"^ 

.ect(A,B) det(X„*r,»*r) 



where we use the notation of Theorem ^.c, and id G Sk is the identity permutation. 
Consequently, 

wini) ■ ■ ■w{Trk)det{K.^) 

(4.6) detiXA.B) = --^^^^^ . 

dei{KintT,intT) W K{a,a) 

aeA 

It is not hard to show H, Lemma 3.1] that for any connected (not necessarily 
planar) resistor network and any subset of vertices V C V, the determinant of the 
corresponding principal submatrix Ky' — {K{a, b))a,b£V' of the Kirchhoff matrix 
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K is positive. Thus all determinants appearing in the right-hand sides of (4.2)-(4.6) 
are positive, and we obtain the following corollary. 

Corollary 4.5. Under the assumptions of Corollary ^.4, (— 1)'^ det(A^.s) > and 
AciiXA.B) > 0. 



The inequality (—!)'"' det(Ayi.s) > in Corollary 45 was first obtained by E. Cur- 
tis, E. Mooers, and J. Morrow Q (for "well-connected circular planar networks" ) 
and by Y. Colin de Verdiere jij (general case) . A short proof based on Theorem LS 
was given in Q 



Applying the inequalities in Corollary 4.5 to arbitrary subsets A' G A and B' C 



B of the same cardinality, we arrive at the following result. 



Corollary 4.6. Under the assumptions of Corollary the submatrices ^Aa.b 
and Xa,b of the hitting matrix X and the negated response matrix —A are totally 
nonnegative. Furthermore, these submatrices are totally positive (i.e., all their 
minors are positive) provided the set Cid{A,B) is not empty, i.e., there exists at 
least one family of vertex- disjoint paths connecting A and B through the interior 

ofn. 

The original motivation for this paper was to provide a combinatorial explana- 
tion for this total positivity phenomenon, in the spirit of Corollary |3.2| . Such an 
explanation is given in Section 0, where total nonnegativity of hitting matrices of 
directed planar networks is established by purely combinatorial means. 

5. Loop-erased walks 

The concept of loop-erased walk (on an undirected graph) has been extensively 
used in the study of random walks, following the work of G. Lawler (see jlj, Sec- 
tion 7] and references therein). Here is the definition, adapted, in the most straight- 
forward way, to the case of oriented graphs. 

Definition 5.1. Let tt be a walk given by 

(5.1) Oq ^ ai ^ 02 ^ ■ ■ ■ ^ a„i . 

The loop-erased part of tt, denoted LE(7r), is defined recursively as follows. If tt 
does not have self-intersections (i.e., all vertices Oi are distinct), then LE(7r) = tt. 
Otherwise, set LE(7r) = LE(7r'), where tt' is obtained from tt by removing the first 
loop it makes; more precisely, find a.; = aj with the smallest value of j, and remove 
the segment — > a^+i-^ ■ ■ ■ ^ Oj trom tt to obtam vr . 

The loop erasure operator LE maps arbitrary walks to self-avoiding walks, i.e., 
walks without self-intersections. In the case of a directed network F = {V,E,w 
this operator projects the weight function on the set of walks, defined in Section 
by w{tt) — Y[ w{ei), onto the corresponding function, denoted by w, on self-avoiding 
walks k: 



w{>c) = w{tt) . 



LE(7r) = > 



Although it will not be needed in the sequel, the author cannot resist a temp- 
tation to mention here a remarkable theorem of Lawler's (first stated in complete 
generality by R. Pemantle |l^). This theorem asserts, speaking loosely, that in 
the special case of a simple random walk on an undirected graph, the probability 
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measure induced by loop erasure on self-avoiding walks from a vertex a to a vertex b 
coincides with the measure obtained by choosing, uniformly at random, a spanning 
tree of the underlying graph, and then selecting the unique path in the tree that 
connects a and b. See for further details. 

6. Minors of walk matrices 

I n th is section, we prove the master theorem on minors of walk matrices (Theo- 
rem state its si mp lified version applicable to walk matrices of planar directed 



networks (Theorem 3^), and give a number of applications. 

The notation and terminology introduced in Section ^ is kept throughout. Thus 
r = (y, E, w) is a directed network with the walk matrix W — {W{a, b))a,bev ■ 

Similarly to Sections || and U, let us choose a pair of totally ordered subsets 
A = {oi, . . . , ak} C V and B — {bi, . . . , bk} C V, not necessarily disjoint, of the 

same cardinality k. Let us denote by Wa,b — {W{ai, 6j))i=i,...,fc the corresponding 

j=i,...',k 

k X k submatrix of the walk matrix W. 

From the definition of the determinant, we have 

(6.1) det(T4^A^i3) = ^ sgn((T) ^ U!(7ri) • • • ■u;(7rfe) , 

"1 > 0(7(1), ■ • • , 0,k ' Ocr(fc) 

where the first sum is over all permutations a &Sk, and the second sum runs over 
all families of walks tti , . . . , TTfe , which connect elements of A to the elements of B 
assigned to them by permutation a. 

Theorem 6.1. The minors of the walk matrix W are given by the formula 

(6.2) detiWA.B) = Y ^gn(^) E 

w(7ri) • • • wiiTk) , 

c<^Sk i < j n LE(7ri) = 



obtained by restricting the second summation in (6.1) to the families of walks 
TTi,...,TTk satisfying the following condition: for any 1 < i < j < k, the walk 
TTj has no common vertices with the loop-erased part of TTi . 

Before embarking on the proof of this theorem, let us pause for a couple of 
comments. 



First, note that the left-hand side of (3.2) is invariant, up to a sign, under changes 



of total orderings of the sets A and B. It follows that the right-hand side possesses 
the same kind of invariance, which is not at all obvious, since the condition imposed 
there on families of walks involves these orderin gs in a nontrivial way. 



Another symmetry of the left-hand side of (6^) that does not manifest itself 
on the right-hand side is the invariance with respect to "time reversal." In other 
words, we can redirect all the edges backwards, while keeping their weights. The 
walk matrix is then transposed, so its minors remain the same, with the roles of A 
and B interchanged. This transformation of the network will replace ordinary loop 
erasure by "backwards loop erasure" (tracing a walk backwards, and erasing loops 
as they appear), and it is not immediately clear that this modification will leave 



expression on the right-hand side of (6.2) invariant (but it will) 



Proof. We will prove (^.2|) by constructing a sign-reversing involution on the 



set of summands appearing on the right-hand side of (6.1), for which the condition 
i < j TTj n LE(7ri) = is violated. The argument will be first presented for the 
special case k = 2, and then extended to a general k. 
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The following general construction will be useful in the proof. Let a — > b be 
an arbitrary walk, and let w be a vertex lying on its loop-erased part LE(7r). The 
edge sequence of LE(7r) is canonically a subsequence of the edge sequence of tt; in 
the latter sequence, let us identify the unique entry of the form v' v which is 
contained in the loop-erased subsequence. (The case w = a is an exception, to be 
kept in mind.) The walk tt is now partitioned at the end of the entry e into two 

walks a v — ^-^ b such that 

(a) the last entry in the edge sequence of tt'{v) contributes to LE(7r'(w)); 

(b) Tr"{v) does not visit any vertices which lie on LE(7r'(u)), except for v. 

The conditions (a)-(b) uniquely determine the partition {n' {v) , n" {v j) of the walk tt; 
in the special case v = a,we set Tr'{v) to be the trivial path v ^ v, and let tt"{v) = tt. 

Let us get back to the proof. Assume k = 2, and let the walks oi 5^(1) and 
a2 ^o-(2) be such that 7:2 and LE(7ri) pass through at least one common vertex. 
Among all such vertices, choose the one (call it v) which is closest to ai along the 

self-avoiding walk LE(7ri). We then partition tti into oi v -^-^ ^(t(i) i following 
the rules above (see (a)-(b)). Let us denote L = LE(7r5^(w)). With this notation, 
condition (b) can be restated as follows: 

(bl) 7r"(w) does not visit any vertices which lie on L, except for v. 

Let us now split TT2 at the point of its first visit to v. More formally, we define the 

partition 02 — ^ v — ^ ^o-(2) of tt2 by requiring that iri^ does visit v before arriving 
at its endpoint. By the choice of v, 

(b2) TTj does not visit any vertices which lie on L, except for v; 

(c) TTj does not visit any vertices which lie on L, except for ending at v. 

Everything is now ready for the path- switching argument. Let us create new 
walks 

1t[{v} TTj 

TTi : fli > V > 6^(2) 

and 

TT2 ■■ 02 > V > 60.(1) ■ 

The map (7ri,7r2) 1— > (7^1,112) is the desired sign-reversing involution. The basic 
reason for this is the similarity of conditions (bl) and (b2), which allows to inter- 
change the portions 7r"(i;) and n!^. It furthermore ensures that the new walk tti 
splits into tt'i{v) = t^i{v) and 7r"(w) = t^2- In particular, the loop-erased part of the 
initial segment of tti remains invariant: LE(7ri(t;)) = L. Now (bl) and (c) show 
that TT2 splits at v, as needed. As a result, applying the same procedure to {tti,tt2) 
recovers (7ri,7r2), and we are done. 

The case of an arbitrary k is proved by the same argument combined with a 
careful choice of the pair of paths to which it is applied. Take a term on the right- 
hand side of ( |6.l| ) which corresponds to a /c-tuple of walks (tti, . . . , tt^). Assume that 
this term does not appear in ( |6.2| ). Thus the set of triples (i,j, w), 1 < ? < j < fc, 
V E V, such that the walks ttj and LE(7ri) pass through v is not empty. Among 
all such triples, let us choose lexicographically minimal, in the following order of 
priority: 

• take the smallest possible value of i; 
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• for this i, let v to be as close as possible to along LE(7ri), among all 
intersections with tTj , for all j > i; 

• for these i and find the smallest j > i such that iTj hits v. 

We then proceed exactly as before, working with the pair of walks (tt^, ttj). Thus 

we partition the walk into — > v - — > 6^(1) i as prescribed by (a)-(b) above. 
Then locate 

• the first visit of vTj to v, starting at aj 

and split nj accordingly into aj — > v — *■ 6cr(j)- Exchanging the portions 7r"(t;) 
and TTj of the walks tt^ and tt^ provides the desired sign-reversing involution. A 
straightforward verification is omitted. □ 



Example 6.2. Consider the directed network in Example 2.1. Let A — {b,c} and 
B = {a, 6}. Then 



3.3) det(W^^,s) = (1 - <Zi<Z2g3)" 



-93(1 - 919293)" 



On the other hand, the only way a pair of walks 7ri,7r2 connecting A and B can 
satisfy the conditions of Theorem 3.1 is the following: take any closed walk b — ^ b 



(so the loop-erased part of tti will be the trivial walk h — >h of length 0), and set 
7r2 to be the walk c — > a of length 1. The corresponding permutation a will be a 
transposition, so sgn((T) = — 1. The weight of 1^2 will be equal 53, and the generating 
function for the weight of vri will be (1 — 919293)^^ • So the right-hand side of ( |6.2D 
will be equal to —93(1 — 919293)^^, in compliance with (|6.3|). 



Example 6.3 (Bernoulli random walk). Let us now look at Example 2.3. The 
results w c arc going to obtain are not new, but the proofs illustrate quite well how 
Theorem 3T works in one-dimensional applications. 

First, let A — {01,02} and B = {61,62} be such that oi < 61 < 02 < 62 and 
02 — 61 = fc , as shown below: 



k 



61 



02 



62 



By formula ( p.Sj ), the left-hand side of (6.2) is 

(p-9)"^ ip-qy^ 



det(W^A,B) 



(p-q) ^q^p ^ (p-q) ^ 



(p-(/)-2(l_qV'). 



Let us now look at the right-hand side. The only possible a is the identity permuta- 
tion. (It might seem that a mistake is being made here, since the scenario oi 62, 
02 ^ 61 may occur without direct collision: the two "particles" can easily cross 
each other's ways without ever being at the same point at the same time. However, 
we only need to make sure that the trajectories intersect — in space, but not neces- 
sarily in time! This distinguishes the loop-erased switching from the conventional 
Karlin-McGregor argument.) 
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Any walk ai bi has the same loop-erased part, namely the shortest path 
from fli to bi . To avoid intersecting that path, a walk 02 62 must never hit 61 . 
The sum of weights of all such walks can be interpreted as 

expected number of visits that a random walk originating 



E{a2, b2,bi) y makes to 62 before hitting 61 for the first time 

(if a trajectory never hits 61, then all visits to ^2 count toward E{a2, 62; ^i))- Thus 
detiWA,B) = [p-qy^il - q'^p-'') = M/(ai, 61)^(02, 62; 61) = (p-(7)"'^(a2, 62; ^i), 

and we conclude that i?(a2, 62; &i) = (p ^ 9)^^(1 ^ q^p^'^) (assuming 5i < 02 < 62 
and 02 — 61 = fc). 

Continuing with this example, let us now swap 02 and 62, as shown below: 



k I 



Ol 61 62 "^2 

Analogous computations give i?(a2,62;fei) — {p — q)~^q^P~^{^ — q^p~^) for 61 < 
62 < 0-2 and 62 — ^1 = ^, 02 — 62 = ^- Finally, let us add another pair of points, and 
change the labelling once again: 

" k I -—> m K 



ai 



hi 



as 



a2 



h2 



Then 



det(w^^,B) = [p-qY 



= (p-g)-3(l-(p'=)(l-(fr). 



6.1 



By Theorem 

(6.4) det{WA,B) = W{ai,bi)E{a2,b2;bi)E{a3, 63; fei, 02) , 

where E(a2, 62; ^i) has the same meaning as before, and 



£^(03, 63; 61, 02) 



expected number of visits that a random walk which 
originates at 03 makes to 63 before hitting either bi or 02 



Substituting W{ai,bi) (p-q)-^ and £^(02, 62; ^'i) = (p- - (f hito 
(|6^), we conclude that 

E{a3, 63; 61, a2) 



(i-(f)')(i-(fr) 



ip-q){l-ii) 



Most interesting applications of Theorem 6.1 arise when the network T is embed- 
ded into a planar domain 51, and the points ai, . . . , and &i, . . . , 6^ are chosen on 
the topological boundary of Q in such a way that the only allowable permutation 
(T is the identity permutation; cf. Figures f|and|4[ 

The following corollary of Theorem 6.1 is useful in such situations. 
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Figure 4. Planar network 



Theorem 6.4. Assume that the vertices oi, . . . , and bi, . . . ,bk are chosen so 
that any walk from a-i to bj intersects any walk from , i' > i, to bji , j' < j. Then 
the corresponding minor of the walk matrix is given by 

(6.5) det(W^A.s) = ^w(^i)---w(7rfe), 

where the sum runs over all families of walks ai — ^- bi, . . . , Ok — ^ bk such that for 
any l<i<k, the walk TTi+i has no common vertices with the loop-erased part ofiTi. 
In particular, if all edge weights are nonnegative, then the matrix Wa,b is totally 
nonnegative. 



Proof. Consider the right-hand side of (3.2). First we note that the loop-erased 



walks LE(7ri) are pairwise vertex-disjoint, and therefore a must be the identity 



permutation, under the assumptions of Theorem 6.4. 

It remains to show that the restrictions on the tt^ may be relaxed as indicated. 
Assume that the walks ai ^- 6i, . . . , Ofe b^ satisfy the condition in Theo- 
rem 6.4. Suppose that the corresponding condition in Theorem 3.1 is violated, that 
is, some walk tt^ intersects LE(7ri), for some i < j — 2. We may assume that among 
all such violations, this one has the smallest value of j — i, which in particular 
means that tTj-i does not intersect LE(7ri). We can then combine segments of TTj 

and LE(7ri) to create a walk a.i bj . By assumptions of Theorem |6.4| , nij must 

intersect the walk flj-i ^'^ . On the other hand, neither LE(7ri) nor TTj 

intersects LE(7rj_i), a contradiction. □ 



Theorem 6.4 can be applied to the study of two-dimensional random walk. Let 



us briefly discuss two possible venues. 

Example 6.5 (Random walk on a planar strip). Consider a Markov chain on the 
state space V = Z x {0, 1, . . . , N} (see Figure ||a) satisfying the following conditions: 

(a) all allowable transitions are of the form (x. ?;) ^ (.x ± 1, y ± 1); 
(6.6) . . ... . '. . 

(b) the transition probabilities are translation-invariant along the x axis. 

Let dTo = : z S Z} and dTi = {{j,N) : j e Z} denote the lower and 

the upper boundaries of the strip. Choose two points a = {i, 0) G OTq and b = 
{i + k,N) £ dTi. Let us denote by T{k) = W{a, b) the expected number of times 
the process passes through 6 provided it started at a. (To simplify matters, let us 
assume that E(k) < -l-oo.) Then the infinite Toepliz matrix T = (T(j — i))ij^z 
is the same as the submatrix {W{a,b))a£dro,bedri of the walk matrix W. Now 



Theorem 6.4 shows that T must be totally nonnegative. Furthermore, it provides 
combinatorial (or probabilistic) interpretations for the minors of T. For example, 
for any positive integers I and m, the 2x2 determinant T{k)T{k + rn — I) — T{k + 
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m)T{k — I) is equal to the expected number of times a trajectory originating at 
(0, 1) hits the point {k + m, N) before hitting {k, N) for the first time. 

The same conchisions will of course hold if the transition probabilities/ weights 
are not translation-invariant, except that the matrix in question will generally not 
be Toepliz. It is natural to ask whether such matrices T can be obtained using 
some modification of the construction presented in Section ^ At present, we do 
not have a good answer to this question. 



hi 



bk 



dTo 



0-3 "^2 0,1 hi 62 &3 



(a) (b) 

Figure 5. Two-dimensional random walks 



Example 6.6 (Random walk on a half-plane). Let us consider random walk on 
the half-plane Z x {0, 1, ... , N} (see Figure ||b). As before, we assume conditions 
(U), and denote (SFo = {(i, 0) : i G Z}. 

Choose two points a — (i,0) € dV^ and b — {i + k,0) G 9Fi, and let H{k) — 
W{a, b). Then the infinite Hankel mat rix T = {T{i+j))i,j^z will be totally nonneg- 
ative, for reasons similar to Example 6.5. Furthermore, it will be totally positive 
provided all transition probabilities between neighbouring states are positive. De- 
tails are left to the reader. 



7. Minors of hitting matrices 



The main goal of this section is to produce analogues of Theorems |6.1| and 
for hitting matrices. As before in Section H, we assume that the vertex set V is 
arbitrarily partitioned into two disjoint subsets: the boundary dT and the interior 
intr. Choose a pair of totally ordered subsets A — {ai, . ..,ak} C V and B — 
C dT, and denote by Xj^g — (X{ai,bj))i^i j; the corresponding 

j=C.'.',k 

hitting submatrix. 

Recall that LE(7r) denotes the loop-erased part of a walk it. 

Theorem 7.1. The minors of the hitting matrix are given by 
(7.1) detiXA,B) - J2 E 

w(7ri) ■ • ■ wiTTk) , 



■^1 I, 

ffli — > K(i 



(1)' 



, Ufe > b, 



'<T(fc) 



where the first sum is over all permutations a ^Sk, and the second sum runs over 
all families of walks tti , . . . , tt^ satisfying the following conditions: 

• TTi has nonzero length, begins at ai, ends at 6cr(i); o,nd in the meantime does 
not pass through any boundary vertices; 

• for any 1 < i < j < k, the walks TVj and LE(7ri) have no common vertices in 
the interior ofT. 



16 



SERGEY FOMIN 



Note that Theorem 7.1 reduces to Theorem 3.1 if the network is acyclic 



Comparing to formula (4.4) (which has a narrower domain of applicabihty) , 
Theorem 7.1 has the advantage of "polynomiality:" the right-hand side of (7.1), 
unlike that of (4.4), is manifestly a formal power series in the edge weights. This 
feature will be essential while extending the result to the continuous case. 

Proof. This theorem can be proved by a direct argument similar to the one 
used in the proof of Theorem 6.1. To save an effort, we will instead use a simple 
observation that will reduce Theorem to Theorem 6.1 



Let us define a new network V by splitting every boundary vertex b e dT into 
a source b' and a sink 6", converting all outgoing edges b—>^c into b'^c, redirecting 
all incoming edges a^b into a^b" , and keeping the edge weights intact. Let dV 
and dV" be the sets of sources and sinks, respectively. Then the hitting matrix 
X of the original network F becomes a submatrix of the walk matrix W' of the 
transformed network F'; more precisely, X{a,b) — W'{a',b") for a, 6 G dV, while 
X{a,b) = W'{a, b") for a €V \dT and b € dT. To complete the proof, it remains 
to carefully reformulate the statement of Theorem 6.1 for F' (with each G dT 
replaced by a- and each bj replaced by b'') in terms of F. □ 



The following analogue of Theorem 6.4 is a corollary of Theorem 7.1 



Theorem 7.2. Assume that the vertices ai, . . . , flfe e V and bi, . . . ,bk £ dT have 
the property that any walk from ai to bj through the interior of T intersects any 
such walk from ai' , i' > i, to bj/ , j' < j , at a point in the interior. Then 



(7.2) 



det(X^^s) = ^ w(vri) • ■ • w{Trk) 



where the sum runs over all families of walks tti , . . . , tt^ satisfying the following 
conditions: 

• TTi has nonzero length, begins at ai, ends at bi, and in the meantime does not 
pass through any boundary vertices; 

• for any \<i<k, the walk 7ri-|_i has no common vertices with the loop-erased 
part of TTi in the interior ofT. 

In particular, if the edge weights are nonnegative, then the matrix Xa.b is totally 
nonnegative. 



The proof of Theorem 7.2 is similar to that of Theorem 3.4, and is omitted 



Example 7.3. In Example uA, let us verify Theorem 7.2. We have 



det{XA,B) = (1 - qiq2q3y 



qiq2qbq& 



qiq^q^qi 
qiq&qi 



= (1 - qiq2q-i) ^qiqiqsqeqj ■ 



On the other hand, the right-hand side of (7.2) is the product of X{ai,bi) = (1 



91 9293) ^9495 and the weight of the only walk 02 62 that does not intersect the 
only self-avoiding walk oi — '-^ 61. This weight being equal to qiqeqr, Theorem 7.2 
checks. 

Hitting matrices of Markov chains. If the weights of a directed network are 



transition probabilities (cf. Example 2.5), Theorem 7.2 has the following probabilis- 
tic interpretation. 

As before, consider a Markov chain on the state space V with "boundary" dT C 
V. The entries of the hitting matrix X are the hitting probabilities: X{a,b) is the 
probability that b is the first boundary state hit by the process that originates at a. 
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(Here and below in Theorem 7.4, a "hit" must occur after the clock is started. In 
other words, if the process originates at a boundary state, it is not presumed to hit 
the boundary right away.) 

Theorem 7.4. Suppose that totally ordered subsets A = {ai,...,afe} C V and 
B = {&!,..., 6fe} C dV are such that any possible trajectory connecting Oi to bj 
through intT intersects any trajectory connecting Oii , i' > i, to bji , j' < j, at 
a point in intT. Then the minor det{XA,B) of the hitting matrix X is egual to 
the probability that k independent trajectories tti, . . . , tt^ originating at oi, . . . , a/j, 
respectively, will hit the boundary dT for the first time at the points bi, . . . ,bk, 
respectively, and furthermore, the trajectory tt^+i will have no common vertices 
with the loop- erased part of -Ki in the interior ofT, for i — I, . . . , k — 1. 

In particular, det{Xji B) > 0, and moreover the matrix Xa b is totally nonneg- 
ative. 

To illustrate, the hitting matrices for two-dimensional random walks discussed 



in Examples 3.5 and |6.6| are totally nonnegative. 



As an application of Theorem 7.4, we obtain another proof of total nonnegativity 



of response matrices of resistor networks (see Corollary 4.6) 



Generalized hitting matrices. For a ^ V and B C dT, let X{a,B) denote 
the corresponding hitting probability; that is, X{a,B) is the probability that the 
process originating at a will first hit the boundary at some point b G B. Using mul- 



tilinearity of the determinant, we obtain the following modification of Theorem 7.4, 
(Other theorems above can also be given similar analogues.) 

Theorem 7.5. Assume that distinct states ai,...,ak G V and disjoint subsets 
Bi, . . . , Bk C dT are such that any possible trajectory connecting Oi to Bj intersects 
any trajectory connecting Oii , i' > i, to Bji, j' < j, at a point in intT. Then the 
determinant det{X (oi, Bj)) is nonnegative, and moreover the matrix {X{ai, Bj)) is 
totally nonnegative. 

Moreover, det{X (oi, Bj)) is equal to the probability that k independent trajecto- 
ries TTi, . . . , TTfc originating at ai, . . . , Ok, respectively, will hit the boundary dT for 
the first time at points which belong to Bi, . . . , Bk, respectively, and furthermore, 
the trajectory Hi^i will have no common vertices with the loop-erased part of tt^ in 
the interior of T , for i — \, . . . ,k — 1. 



In order for the statement of Theorem 7.5 to make sense, the Mark ov p rocess 
under consideration does not have to be discrete. For example, Theorem [7.5| gener- 
alizes straightforwardly to (non-isotropic) Brownian motions on planar domains, or 
arbitrary simply connected Riemann manifolds with boundary. (See Figure H.) The 
proofs can be obtained by passing to a limit in a discrete approximation. The same 
limiting procedure can be used to justify well-definedness of the quantities involved; 
notice that in order to define a continuous analogue of the probability appearing in 



Theorem 7.5, we do not need the notion of loop-erased Brownian motion. Instead, 
we discretize the model, compute the probability, and then pass to a limit. One can 
further extend these results to densities of the corresponding hitting distributions. 
Technical details are omitted. 

The rest of this section is devoted to a couple of characteristic applications 
involving Brownian motion. 



18 



SERGEY FOMIN 



ai 03 
02 



Figure 6. Theorem 7.5 



Brownian motion in the quadrant with reflecting side. Let 

n^{ix,y) : x>0,y>0}, 

and consider the Brownian motion in Q which reflects in (bounces off) [the positive 
part of] the x axis (see Figm-e|7^). The y half-axis is absorbing, i.e., the process 
stops once it reaches a state {0,y), y >0. This is equivalent to taking a Brownian 
motion in the half-plane a; > that stops upon hitting the y axis, and reflecting 
the portions of it that go into the lower quadrant {{x,y) : a; > 0, y < 0} over the 
X axis. 



yi« 



yi y2 y 

(I 

X X 



Xl Xl Xl Xl 

(a) (b) 

Figure 7. Brownian motion in domains with reflecting sides 



Let the process start at time t = at a point (a;,0), a; > on the x axis. It 
is classically known (see, e.g., Section 1.9] for two different proofs) that the 
Brownian motion in half-plane hits the y axis for the first time at a point Y[x) that 
has the Cauchy distribution with appropriate parameter. Taking absolute values, 
we conclude that our process in O induces the hitting distribution on the half-line 
{(0,2/) : y> 0} which has the density 

1x 

K{x,y) = 12, — 2T' y^O- 
■n{x^ -f y^) 

It follows from the general results above that the kernel K{x, y) is totally positive, 
i.e., any determinant of the form (lei{K{xi,yj)) with Q < xi < X2 < ■ ■ ■ and 
< yi < 2/2 < ■ • • is positive. This result is not new. What seems to be new is 
the interpretation we give to such determinants in terms of loop-erased walks. To 
illustrate, let us fix < a;i < X2-, and consider two independent trajectories tti and 
1T2 which start at points (a;i, 0) and (a;2, 0) and end up hitting the y axis at points 
(0, Y{xi)) and (0, Y{x2)): respectively. Assume that these two points are the points 
(0, yi) and (0, 2/2) with < 2/1 < 2/2, but we do not know which trajectory hit which 
of the two points. Then the conditional probability of the scenario Y{xi) — yi, 
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Y{x2) = 2/2 is larger than that of Y{xi) = 2/2, Y{x2 
these conditional probabilities, which is given by 

K{xi,yi)K{x2,y2) - K{xi,y2)K{x2,yi 



yi, and the difference of 



K{xi,yi)K{x2,y2) + K{xi,y2)K{x2,yi) ' 

is equal to the conditional probability that 1:2 did not hit the loop-erased part of ni . 
(We repeat once again that the well-definedness of the loop-erased Brownian motion 
does not have to be justified in order for such probabilities to make perfect sense.) 
Computations show that 



K{xi,yi) K{xi,y2) 
K{x2,yi) K{x2,y2) 



AxiX2{xl - xDiyl - yl) 



and therefore the conditional probability in question is given by the surprisingly 
simple formula 

{xl - xl){yl - yl) 



(7.3) P(^i nLE(7r2) = I {Y{x^),Y{x2)]^{yi,y2}) = 



If we do not condition on the locations of the hitting points (still keeping the 
initial locations xi < X2 fixed), then the probability that the trajectory 7:2 origi- 
nating at X2 will not intersect the loop-erased part of the trajectory tti that starts 
at Xl (or similar quantity with tti and 7r2 interchanged) is equal to 



P{Y{xi)<Y{x2))-P{Y{xi)>Y{x2)) = 



VI 



K{xi,yi) K{xi,y2) 
K{x2,yi) K{x2,y2) 



Computing the integral yields 
P(7ri nLE(7r2) = 0) = 



Li2 (-a) +U2{l~a)+ ln(a) ln(l -I- a) -|- 



dy2dyi 



12 / ' 



X2 

where a — — and Li2(t) =7 ^ is the dilogarithm function. In particular, if 

n—l 



_ \/5+l 



(so the initial locations xi and X2 are in golden ratio), then we obtain. 



using [|, the formula P{ni n LE(7r2) = 0) = i - -^(ln(^^))l 

Brownian motion in a strip with reflecting side. In this example, we let 

O = {{x,y) : a; e R,0 < y < 1} 

and consider the Brownian motion in which is reflecting in the x axis (see Fig- 
ure 0b). The process begins at a point {xq, 0) on the x axis and stops when it hits 
the line y = I. 

Standard computations (employing either conformal invariance of Brownian mo- 
tion or the reflection principle) yield the well-known formula for the hitting density: 



e 2 



2cosh(f^) 



Thus A" is a totally positive translation-invariant kernel; equivalently, the function 
K{x) = ^ is a Polya frequency function. This particular example of a 

PF function is of course well known; see, e.g., Karlin ||l^, § 7.1, example (b)]. The 
computation of probabilities related to loop-erased walks in this model is omitted. 
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